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^ ! Abstract 

Full quantum dynamics calculations have been carried out for the ionic reaction 4 He^ + 3 He and 
' state-to-state reactive probabilities have been obtained using both a time-dependent (TD) and a 

time-independent (TI) approach. An accurate ab-initio potential energy surface has been employed 
for the present quantum dynamics and the two sets of results are shown to be in agreement with each 
other. The results for zero total angular momentum suggest a marked presence of atom exchange 
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(isotopic replacement) reaction with probabilities as high as 60%. The reaction probabilities are 
only weakly dependent on the initial vibrational state of the reactants while they are slightly more 
sensitive to the degree of rotational excitation. A brief discussion of the results for selected higher 
total angular momentum values is also presented, while the Z-shifting approximation 1] has been 
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used to provide estimates of the total reaction rates for the title process. Such rates are found to 



be large enough to possibly become experimentally accessible. 
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I. INTRODUCTION 



The study of helium nanodroplets has been shown over recent years to provide a very 
interesting and novel medium that acts as a "quantum matrix" in which one can probe the 
ro- vibrational and electronic spectroscopy of the various dopants which can be "solvated" in 
it and in which one can also observe the nonclassical effects caused by the surrounding He 
atoms , together with the changes due to their bosonic/fermionic symmetry properties 
^J. The experiments which can perform ionization of the helium matrix further raise the 
issue of the mechanism leading to ionization vs fluorescence, and of the understanding of 
the features of the interaction of the resulting ion with the atoms in the droplet which can 
show evidence of marked quantum many-body effects i| . 

One of the nanoscopic consequences of the formation of a permanent cationic impurity 
within the droplet is its high mobility within the weakly interacting solvent, whereby the 
positive charge can in principle migrate within it by some resonant charge hopping mech- 
anism []|. This process is taken to be terminated either by formation of a strongly bound 
let ion or by charge transfer to a dopant species. Our recent calculations on that process 
have shown that termination can also be achieved by multiple inelastic collisions within 
the droplet that enhance the evaporation of atoms from the latter. Furthermore, possible 
branching reactions of either the dopant ions or of the ionic moiety He^~ could also occur 
within the droplet because either species will trend to move towards the center of the cluster 
in order to minimize the total potential energy of the composite system jljj. 

Recent experiments have been carried out on specific reactions initiated by primary ion- 
ization of the droplet 12 , 3| and have shown that a rapid quenching of the internal degrees 



of freedom of the partner species is favored by the superfluid environment. The correspond- 
ing reactive branching of the ionic partners is therefore also markedly affected and could 
cause, in some cases, stabilization of intermediate complexes different from those expected, 
and seen, in the gas phase 3]. In a recent theoretical and computational study [14] we have 
in fact shown the possibility of having strongly exothermic chemical reactions without acti- 
vation barrier, the occurrence of which should be greatly facilitated by the presence of the 
inert, quantum environment at the low temperatures existing within the helium "nanocryo- 



stat" 



In the present analysis we therefore wish to present calculations where the stabilized 
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FIG. 1: Coordinates used in the TD code 

primary ion in the droplet, the He^~ species occupying its lowest ro- vibrational level, is made 
to react with one of the adatoms of the environment. The idea is to analyze several aspects 
of the problem which carry considerable theoretical interest and which can also guide us to 
understand the possible behavior of the ionized droplet under actual experimental conditions. 

In particular, taking advantage of our recent calculations of the potential energy surface 
(PES) for the He^+He system [9|, we have studied the quantum reaction for the isotopic 
exchange, i.e. for the molecule 4 He^~ colliding with neutral 3 He, in order to add physical 
distinguishability within the process at hand. We have carried out the calculations, as we 
shall further describe below, using both a quantum time independent (TI) approach and 
a quantum, time dependent (TD) wavepacket approach. As we shall see, both methods 
produce essentially the same results and therefore provide a sort of internal check for our 
findings thereby adding further credibility to the present theoretical "experiments". 

II. THE TD CALCULATIONS 

In the 3 He + 4 He^(t>,j) collisional event two identical bosons 4 He(S=0) forming 
a diatomic ion collide with a different isotope (fermion) 3 He giving rise to inelastic 
(— > 3 He+ 4 He^(t> ', j')) and reactive (— > 4 He + 3 He 4 He + (t>', j')) processes. 

In a Body-Fixed (BF) frame of reference defined with respect to the Laboratory frame by 
the three Euler angles 9, <p, x, we have chosen to use bond coordinates as the three internal 
coordinates describing the relative position of the three particles (see Figure 1) that is two 
atomic distances R\, R2 and the angle a between them. A method based on these coor- 
dinates was presented previously and applied to the calculation of state-to-state reaction 
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probabilities because they would require relatively complicated changes of coordinates. Our 
present choice of bond coordinates allows a relatively simple and fast evaluation of these 
state-resolved probabilities when only two molecular arrangements are accessible of the three 
which exist in a generic ABC system. The method is applied here for the first time to a 
case in which all the three arrangements are open but where we can still exploit the AB 2 
symmetry for the calculation of state-to-state probabilities, as we shall further explain later 
in this section. 

A convenient BF frame is the one defined such that the z-axis lies along the vector joining 
the center of mass of the initial diatom to the 3 He atom (see Fig. 1) and the three atoms 
are in the xz plane. In the chosen set of coordinates (9, 0, Xi Ri-> Ri-> at), the complete 
Hamiltonian for an arbitrary value of the total angular momentum J takes the form [1 01] 



H = -*L f^-^-R 2 ^) - *1 f^-^-R 2 ^] + f 1 4- 1 ^ f 2 + 
2^ \R 2 8R l 1 dR 1 J 2/i 2 \R 2 0R 2 2 dR 2 J V 2 /^? ^ 2 R 2 J 

+ (J2 2 ~^" J) + f 12 + V(R l , R 2) a) 



;i) 



where we denote by m , mi and m 2 the nuclear masses, with tuq being that of the "reference" 
atom, hence fii = m Q mi/ (rriQ+m^, fi 2 = m m 2 / (m +m 2 ) are the reduced masses associated 
with R,! and R2, while /i = m 2 (m + mi)/(m + m x + m 2 ) is the reduced mass associated 
with the R Jacobi vector. 

In Eq.(£Q) L 2 is the angular momentum operator given by the following expression: 

#2 .2 f 1 d . 8 1 d 2 } 

sin a oa oa sin a ox ) 

while J is the total angular momentum operator. The term ^J 2 — 2J • jj /2fiR 2 where j 
is the angular momentum of the diatomic molecule in the reactants' Jacobi coordinates, 
incorporates the Coriolis couplings given by: 

i • f ( QRi \ 9 d n ■ 9 d 
J ■ j = -in < — — : cot a J x - qRi sin aJ y —— + 

[\R 2 sma J dx oR 2 

( qR\ cosaA ±9 ; 9 1 

where q = fii/m . 

Finally, in Eq. ((TJ T 12 is the kinetic energy cross term due to the non-Jacobian character 
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FIG. 3: Lowest ro-vibrational levels of the 4 He- 4 He (left) and 3 He- 4 He(right) molecules. For the 
former we have reported only the j = 1 levels while for the latter we include only the j = 0. 
Energies are in eV, measured for the bottom of the entrance valley of reaction 

of the internal coordinates we have chosen and has the form 

h 2 { sin a d d 2 sin a d 2 sin a d 2 



T\2 — — { - - » cos a-— — — — H — — — „^ - — h 



m [ R\R 2 da dRidR 2 R 2 dR\da R\ dR 2 da 

cos a 



(4) 
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m RiR 2 

When using the bond coordinates we can distinguish three regions in the PES associated 
with three different arrangements as described in figureEl the entrance arrangement (that we 
shall call the channel 1) [ 4 He 4 He] + + 3 He and the two other possible product arrangements 
(2 and 3) that describe the [ 4 He 3 He] + + 4 He and [ 3 He 4 He] + + 4 He combinations. The 
arrangements 1 and 2 correspond respectively to coordinates R% or R 2 remaining finite 
while arrangement 3 would lie along the direction a = as both R\ and R 2 become very 
large. The inelastic process will be associated with the species remaining in channel 1, while 
both channels 2 and 3 will correspond to the reactive process: although, apparently treated 
in a very different way due to our system of coordinates, these two channels are physically 
equivalent. 



6 



A. Wavepacket representation 

The total wavepacket is, in the general case, expanded as[ic| 

* JM *(9,cf>, X ,R h R 2 ,a,t) = J2 W&iiMx) ^ R ^ a ^ 

ti\ti 2 



(5) 



where W^f n are linear combinations of Wigner rotation matrices|15] of a given parity under 
inversion of all coordinates, e, and where M and Q are the projections of the total angular 
momentum, J, on the space-fixed (SF) and body fixed (BF) z-axis, respectively Insertion 
of Eq.(jHJ) into the time-dependent Schrodinger equation using the Hamiltonian of eq. (PJ), 
generally yields a set of first order differential equations for the $Q Me (i?i, R2,a,t) coefficients. 
For the case of J = which will be considered in the present calculations the previous 
expansion contains only one term. The wavefunction is then given by the solution of the 
equation 



h 2 d 2 



h 2 d 2 



2/xi dR\ 2fi 2 dRl 
with the following meaning of l 2 and t\ 2 



+ 



2^! + 2^) P + t - + l/ }* 



d 



d 



— sin a — 

sm a oa oa 



+ cos a 



tl2 

d 



+ 



1 







sma 
d 



d 1 

+ 







da \R1dR2 R 2 dR! 
1 d 2 



+ 



Ri (9i?2 R2 dRi R1R2 dR\dR 2 



R\R 2 
cos a 
m R 1 R 2 



(6) 



-1- 



The integration of the above equation was performed by using the Chebyshev method 
and the &(Ri, R 2 , i) coefficient is represented on finite grids for the internal coordinates 
Ri, R 2 , a. A set of equidistant points, R\, R J 2 , was chosen for the rectangular bidimensional 
radial grid fni xn 2 ), which allows the evaluation of the radial kinetic terms using Fast Fourier 
Transforms (FFT) |l7j . For the angle a a set of n a Gauss-Legendre quadrature points, a k 
(with weights u k ), is used. Thus, the terms involving derivatives in a are evaluated through 
a discrete variable representation (DVR) transformation which reduces the procedure to a 



simple multiplication of a matrix by a vector 
of the wavepacket is then given by 
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23| . The grid representation 
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TABLE I: Parameters used in the wavepacket propagations for v = 0, 1, 2 j = 1 and J = 
0; distances in A, times in ps. The gaussian parameters (JCq,T) in Eq (9), are given by 
(14.30, 0.42), (10.87,0.55), (6.27,0.96) respectively 
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0.69, 19.50, 360 


( Tyvnin ( °A \ T~>'tn(ix { °a \ \ 

{K 2 ' (A), H 2 ' (A), ri2) 


0.69, 19.50, 360 
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14.10, 0.016 


(Rf s (A),T 2 ) 


14.10, 0.016 




13.00, 13.00 


At(ps) 


0.003 



where, for convenience, the Gauss-Legendre weights are introduced. 

In order to use a finite bidimensional radial grid, the wavepacket is absorbed at each time 
step by multiplying the wavepacket for / 1 (i? 1 )/ 2 (-R2), where fi(Ri) = exp[—Ti(Ri — R°; bs ) 2 ] 
for Ri > R1 bs and fi(Ri) = 1 otherwise. The absorbing regions are presented in the lower 
panel of Figure El shaded area. 

The actual parameters of the propagation used in the calculations are listed in Table 
I. It should be noted that the number of angles, n a , required to converge is indeed very 
large: channel 3 is asimptotically open only for a = and many angular quadrature points 
are required in order to have a grid dense enough in the small-angle region. The initial 
wavepacket represents the reagents approaching the collision region with the diatom in a 
given (z/, j) ro- vibrational state and with a continuum distribution of relative kinetic energies. 
The initial wave packet can thus be expressed (in reactants' Jacobi coordinates (r,R, / y)) as 
a product of a diatomic wavefunction and a Gaussian function for the relative translational 
coordinate: 

Ht = 0) = X „Ar)Y jno (l,0)G(R) (8) 

where G(R) is a complex Gaussian function written as follows 

1/4 r (R-Ro) 2 



r 2 



— iJC (R — Rq 



(9) 
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The gaussian is centered at a convenient Ro value such that the interaction between the 
reactants can be considered negligible. The energy distribution for the initial wavepacket 
a(E) is then determined using 



This system contains two identical bosons and hence the total wavefunction must obey 
the correct spin statistics: the total wavefunction must have the correct symmetry under the 
action of the P\2 operator that exchanges the two bosons. This means that (as in all the AB 2 
systems) the calculations with even and odd initial j values for the diatomic wavefunction 
are decoupled: since the Hamiltonian commutes with the parity operator P12, a wavepacket 
built as an eigenfunction of Pu (as it occurs with our initial wavepacket <3>(t = 0)) will 
remain so under time evolution. An additional simplification in our calculations stems from 
the consideration that, since the 4 He^" is a molecule composed of two spinless nuclei, 
only odd rotational states are allowed in the entrance arrangement. As a consequence of 
these constraints, it is clear that the appearance of the forbidden even-j rotational states in 
the inelastic distributions obtained from the evolved wavepacket using eq. (11) below would 
indicate convergence problems. This was one of the various criteria we used for checking 
convergence in our present calculation and we found no instance where such states would 
appear. 

The main steps of the TD method applied to the J=0 case were the following jlp| |: 

• The initial wavepacket, obtained in the reactants' Jacobi coordinates (r, R, 7) is trans- 
formed into bond coordinates (Ri, R 2 ,a); 

• The propagation is performed using bond coordinates; 

• The S matrix elements for the transition into any other reactive or inelastic final 
state can be obtained from the evolved wavepacket using the product asymptotic 
analysis method of Balint-Kurti et al. (24], 13 that requires a change of coordinates 
to the Jacobi system (R'^, r', 7') corresponding to the final arrangement. From the 
transformed wavepacket at each time step, the state-to-state probabilities are obtained 




(10) 



with , k, 
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using the expression 

where, as noted above, k u if = J 2fi(E — E„ij/)/h 2 , fi is the reduced mass of the 
desired arrangement channel, R'^ is a value of the Jacobi radial coordinate well into 
the asymptotic region and the A^j/^R'^, E) are the energy transforms defined for the 
desired channel as: 

POO 

A u , j ,(R' 00 ,E)= / rfte^(x^(0^o(cos 7 ')|$(^,r', 7 ',t)) (12) 
Jo 

The computational cost of the required coordinate transformations is not very high be- 
cause the arrangement channels 1 and 2 share one of the Jacobi radial coordinates with the 
bond coordinates system r' = Ri, thereby reducing the numerical effort due to coordinates 
tranformations: even in the general case of J ^ 0, only simple rotations on the plane identi- 
fied by the three atoms would be required to move to the desired BF system. Arrangement 
3 would require a much more expensive change of coordinates, but due to the symmetry, 
the probability for the 3 He to strip either of the two 4 He atoms must be exactly the same. 
Hence, we can assume that the probabilities in channel 3 are the same as those in channel 
2 and thus we can avoid calculating them; any deviation from unity of the sum of all of 
them can then be taken as a sign of lack of convergence. This has been another criteria in 
determining the convergence of our calculations in the actual situations of this work. 



B. The TI Calculations 

We have employed the abc program of Skouteris et al. |26( which performs an expansion 
of the total wavefunction in Delves hyperspherical coordinates [22J and uses the coupled 
channel method formulated by Schatz |2J| with an additional orthogonalization scheme of 
the vibrational basis at fixed hyperradius to avoid over-completeness of the basis set in the 
short range region j^J. The coupled channels equations are then solved using the constant 



reference potential log-derivative algorithm of Manolopoulos 30] . 

The convergence of the resulting S matrix is sensitive to two parameters, the hyperradius 
Pmax at which the asymptotic condition is imposed and the step-size Ap of the propagator. 
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Converged results within 1% were obtained with Ap=0.02 a.u. and p max = 25.0 a.u. for the 
lowest collision energy considered here (~7 meV). 

The basis set was built using all the ro- vibrational levels of the two [He-He] + diatomics 
whose energies lie below the cut-off energy value of 1.5 eV as measured from the bottom of the 
entrance channel. Since our calculations would have otherwise involved the use of too many 
rotational levels, we have imposed a maximum value of j ma x=20 for each molecule which 
is enough to obtain probabilities converged within 5% also for the highest rotational states 
included. The basis set chosen has been able to reproduce almost perfectly the TD results 
also at the higest energies employed here as we shall see below. Our basis set comprised of 
a total of 235 basis functions for J = 0. As noticed before the [ 4 He- 4 He] + molecule can only 
have odd rotational states. 

We have performed the calculations for about a thousand different collision energies on a 
grid ranging from 7.375 meV up to 1 eV. The lowest collision energy has been chosen to be 
1 meV higher than the energy necessary to make the reactive channel for the title reaction 
to be energetically accessible when considered in the following direction 

4 He J" + 3 He — > 3 He 4 He + + 4 He (13) 

The small endothermicity of 6.375 meV is due to the difference of zero point energy 
between the v — 0,j — 1 state of the reactants in eq. (j!3)l and the v = 0,j = state 
of the heteronuclear product molecule. For clarity we report in Figure El a diagram of the 
vibrational energy levels involved in the calculations. 

III. RESULTS AND DISCUSSION 

In Figure 0] we report the reaction probability when starting with the initial state 
[y = 0,j = 1) and summing over all the final accessible states. The two sets of data 
that refer to the two different, independent calculations we have mentioned above are in 
good agreement over the whole energy region of interest and provide a useful internal check 
of our calculations. The small differences are due to different causes: first the truncated size 
of our basis set expansion in the TI code is the source of a small level of inaccuracy espe- 
cially for the transitions involving the higher rovibrational states; in second place, the TD 
calculation become inaccurate (as expected in general) for very low kinetic energies (below 
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10 meV in the present case) because of the finite time for progation that necessarily neglects 
the contribution to the correlation functions coming from longer times. The abc program 
uses in fact a basis set which is constructed as an orthogonalized product of ro-vibrational 
arrangement wavefunctions. As a consequence, the description of reactions where triangular 
geometries are dominant (e.g. Li2+Li system) may be inaccurate. Although our system here 
las clearly shown a transition state and a minimum energy path that is mainly collinear 
El, we still intend to provide a consistency check of our calculations by comparing the 
TI results with those from a TD method that is not based on basis sets expansion. On 
the other hand, the TD method is based on a time propagation in bond coordinates and 
might therefore present some problems when applied to a A-B2 system where the third ar- 
rangement (channel 3 in Figure EJ) becomes energetically accessible when a ~ 0; this would 
require many angular points in order to describe the reactive process as discussed earlier on 
in this work. 

The pattern of resonances that can be seen in the reaction probabilities profile is very 
complicated as many of these features may be due to the opening of the various asymptotic 
channels of the two diatomics involved in the reaction. Many others are likely to be due to 
Feshbach-type resonances with metastable states of the triatom and to potential trapping 
due to centrifugal barriers in the exit channels, but it would be premature to analyze them 
now in any detail given the absence of some experimental indication as to their physical 
presence. 

In FigureElwe report the individual contributions of the reaction probabilities for different 
final vibrational channels (the probabilities are rotationally summed). As can be seen from 
that figure, the total reaction probability is larger than 60% over most of the range of collision 
energies considered here (the collision energy in this figure is measured with respect to the 
4 He 2 (^ = 0, j = 1) channel i.e. to the reactants' ground state). However, at large collision 
energies the probability of populating the ground vibrational state of the products decreases 
and the flux is more or less equally redistributed over all the open, final vibrational states. 

In the highly symmetrical systems like A3 one should expect that the total reaction 
probabilities would be exactly the same as the non-reactive one simply because the two 
fluxes are not distinguishable. In our case, however, a certain degree of asymmetry is 
introduced by the isotopic change and therefore the reactive process becomes more likely to 
occur than the inelastic collision process. This feature, at high energies (above 1.0 eV) is 
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FIG. 4: Reaction probabilities for TD (dotted lines) and TI (solid lines) calculations for the initial 
state [y = 0, j = 0). The lower panel shows the same probabilities on a logarithmic scale for the 
lowest energies above reaction threshold 

however mainly due to the fact that the there are two arrangements containing the product 
AB molecule. At these energies, in fact, the helium exchange reaction has a total probability 
of roughly 60%, 30% for each of the two possible exchange, while we find a 40% probability of 
simple elastic/inelastic non-reactive scattering. This behavior may be attributed to the mass 
difference because it seems reasonable that when the collision takes place at high energies 
the lighter 3 He may be not efficient in substituing the heavier isotope. We also believe that, 
given the relatively high collision energies employed here, the small difference of zero point 
energy mentioned in the previous section has a very small effect on the dynamical behavior. 
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Collision Energy (eV) 



FIG. 5: TI (solid lines) and TD (dashed lines) rotationally summed reaction probabilities for 
different final vibrational states. The initial state of the reactant molecule is {y = 0, j = 0). 

Below 1.0 eV we see from Figure H] that the probability may rise well above 0.6 and reaches 
the range 0.8-1. In this case the probability of a non-reactive scattering is less than 20%. 
We believe that the dominant effect at these lower energies may be the higher density of 
states of the product AB molecule with respect to the reactant A 2 molecule which has only 
odd rotational states for the symmetry reasons mentioned before. 

Much more difficult is the task of finding significant patterns of behavior in the final 
products' rotational distributions. The first three state-to-state reaction probabilities for 
different final rotational states, when the reaction begins with the ground state reactants 
are reported in Figure El We can immediately see there that high degree of variation of each 
individual probability persists as j changes as it also does for the other open channels not 
shown in the Figure. At higher energies (not shown in the Figure), on the other hand, the 
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Collision energy (eV) 

FIG. 6: A limited portion of the TI state-to-state reaction probabilities for ground state reactants 
and various final rotational states i.e. for the reaction 4 Hejj~(f = 0, j = 0) + 3 He — * 4 He 3 He + (V = 
0, j') + 4 He where the final molecule is always in the v = vibrational level 

various final rotational states of the product molecule are all substantially populated. 

Another important piece of information that we may obtain from our calculations is the 
effect of increasing the internal energy (i.e. the ro-vibrational excitation) of the reactants. 
In Figure [3 we therefore report the total reaction probability summed over all the open 
final states but now for three different initial states [y = 0, j = 1), [y = l,j = 1) and 
[y = 2, j = 1) (the collision energies have been selected relative to the specified initial state). 
As it can be easily seen for this figure, the initial degree of vibrational excitation (at least 
for the lowest states examined here) does not change significantly the overall dynamical 
behavior of the system: the reaction probabilities rise rapidly at threshold from zero to 
~ 60% and remain roughly of that size when the collision energy increases. 

15 
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FIG. 7: TI reaction probabilities for different initial vibrational states (solid lines); TD results are 
also reported as dashed lines 

A similar situation has been obtained when considering the reactions probabilities pro- 
duced by rotationally excited reactants. For example when looking at the results reported in 
Figure |Hl we see how even a large rotational energy content does not produce a corresponding 
increase in reaction probabilities. It is worth pointing out however that we are limiting this 
analysis to the J = case and therefore a situation in which the initial molecular angu- 
lar momentum is balanced by a corresponding relative orbital angular momentum I. This 
means that the reaction paths that we sampled here with rotationally excited molecules 
would present centrifugal barriers. In conclusion, however, we can say that the reaction 
dynamics in our system is probably dominated by resonances and therefore the memory of 
the initial rotational state is soon lost during the reaction. 
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Collision Energy (eV) 



FIG. 8: TI reaction probabilities for different initial rotational states (the reactant molecule is in 
the v = vibrational state) 

IV. HIGHER ANGULAR MOMENTA 

To increase the value of the total angular momentum beyond J = for reactive scattering 
is still a computationally demanding task. From the analysis of the lower partial waves, 
however, we can already extract useful information about the overall dynamics. When going 
to the situation with J ^ the number of coupled channels increases due to the proliferation 
of the possible VL values, where Vt is the helicity quantum number j26(. Furthermore, for 
each J > value we have two possible values of the total parity eigenvalues (p = ±1) which 
require two independent calculations. The calculations for J ^ that we present here are 
"exact" (no approximations such as coupled states have been used) and have been done 
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using the TI method outlined above and the abc code. 

We report in Figure I§1 and ITU1 the total probabilities for J = 1 and J = 2 for the reaction 
which starts with 4 He^~(V = 0, j — 1). In both cases the results are compared with the J = 
reaction probabilities (thick solid lines in the upper panels). As can be seen in the upper 
panels of those figures the reaction probabilities for higher angular momenta are very similar 
to the ones for J = 0, the only difference being the small shifts in the resonance pattern. This 
is true also at relatively low energies although the resonance pattern is becoming different 
especially close to threshold. 

When looking at the lower panels of Figs. El and El we can further see the contributions 
due to the initial Q — 1 helicity component of the reaction: as it may be expected for a 
reaction which has a marked collinear constraint this last contribution is significantly lower 
than for Q = especially at lower energies. 

Given the fact that accurate calculations for other values of J > would require com- 
putational times which are too long and since we are here in presence of a ionic potential 
which in principle may require many total J in order to properly converge to a reactive cross 
section, we have decided to calculate the reaction rate constants by using something similar 
to a J-shifting approximation, i.e. by using only the J = reaction probability P(E) to 
estimate the total rate constant. A realistic approximation when dealing with barrierless 
(ionic) systems like the one we are examining here, has been suggested in Ref [lj: this ap- 
proximation essentially consists in obtaining the J ^ reaction probabilities by judiciously 
"shifting" the J = values. We have here used the formula: 

Pj(E;v,j)~P J=0 (E-VJ;u,j) (14) 

where u, j identify the initial state of the reactants and Vj is the height of the centrifugal 
barrier of the entrance channel taken along a suitable monodimensional potential generated 
by a non-zero initial orbital angular momentum I value that is allowed for a given J and j; 
we have chosen to use the potential of a collinear geometry given by He — r eq — He — r — He 
with r m = 2.046 a.u.. The only difference with ref. Q is that we have initially a j = 1 
molecule which means that for each J value there are three allowed values of / given by 
| J — 1| < I < J + 1. Since this happens for all J values that may contribute, the further 
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FIG. 9: TI reaction probabilities for J = 1, Q = (upper panel) and for J = 1, = 1 (lower 
panel) 

initial-j averaging leaves the cross section unchanged 

a(E coll ; u, j) = ^- ^(2 J + l)Pj(E; u, j) (15) 

K coll J 

From this cross section we finally obtain the rate constants reported in Figure ^2 Al- 
though this is a very approximate prescription, our final result can still tell us that the rate 
constant for a reactive exchange process in Hejj" is of the order of 10~ n — 10~ 9 cm 3 • s _1 
a value that is significantly lower than the Langevin capture rate also reported in Figure 
im One should note also that, although Langevin capture rates assume the reaction to be 
exoergic, the very small endoergicity of the present system still allows us to use it for an 
estimate of the reaction rates. 
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FIG. 10: TI reaction probabilities for J = 2, 17 = (upper panel) and for J = 2, O = 1 (lower 
panel) 

V. PRESENT CONCLUSIONS 

We have presented new theoretical results on the reactive dynamics of He^ system. In 
particular, we have discussed the reactive behavior for the lowest total angular momenta 
and have extended the calculations to obtain approximate rate constants for the title re- 
action. In order to make the process physically clearer and to get a better insight into its 
mechanism, we have substituted one of the 4 He with an 3 He so that we have introduced 
a small endothermicity in the atom exchange process that is occurring during the reactive 
event. In this way we have also introduced an asymmetry in the density of states associated 
with the reactants (here an homonuclear molecule) and the products (here an heteronuclear 
molecule). This difference makes for a more favorable reactive process with respect to to 
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FIG. 11: Rate constants for the reaction 4 He 2 + 3 He — ► 3 He 4 He + 4 He as a function of temper- 
ature. Also shown as an horizontal line is the Langevin value. 

the simpler inelastic collision process as shown by our calculations. 

The reaction considered here has a collinear MEP without activation barrier and therefore 
behaves in the main like a typical ionic reaction. We have obtained complete state-to-state 
probabilities for the reaction using a time independent method and then, because there may 
be some issues related to the use of an arrangement-based basis set expansion, we have 
checked the numerical reliability of our TI findings by further reproducing our results using 
a time dependent procedure. 

Our present results clearly show that the reaction under study turns out to be quite 
efficient when J = and represents more than 60% of the scattered flux at the energies con- 
sidered. Internal excitation of the colliding partner does not appear to produce substantial 
increases of the reaction probabilities at least for the lowest vibrational and rotational states 
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of the reactants. Even with relatively highly rotationally excited He^ (up to j = 19) the 
size of the reactive probability remains very similar to the one for the non rotating molecule. 
This may be due to the fact that the reaction mecanism is going through many resonances, 
thereby dynamically loosing the effects of having "hot" reactants. 

Furthermore, by looking at the final distributions over the vibrational channels we have 
seen that the reaction can easily produce vibrationally excited molecules and the reactive 
flux is more or less equally redistributed over all the open, final vibrational states. The final 
rotational populations, as far as we can judge from our data, seem to be strongly dependent 
on the collision energy and do not show any simple pattern of interpretation. 

We have also performed fully converged calculations for J = 1 and J = 2 and we have 
seen that the latter data provide very similar results to J = calculation although the 
resonance patterns slightly change and shift in energy, as one expects from such systems. 
We have thus been able to provide rate constants for the reaction through the use of a 
/-shifting approximation which yields final rates of the order of 10~ 10 cm 3 • s^ 1 . The 
latter results indicate that in the present system the reactive process may be as efficient 
as the inelastic vibrational de-excitation process estimated earlier by us p| in providing a 
mechanism of energy exchange in the droplets. The size of the exchange rate also suggest 
that observation of such reaction may be possible in the droplet environment. 
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